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ABSTRACT 

We investigate a stationary pair production cascade in the outer magnetosphere of 
a spinning neutron star. The charge depletion due to a global current, causes a large 
electric field along the magnetic field lines. Migratory electrons and/or positrons are 
accelerated by this field to radiate curvature gamma-rays, some of which collide with 
the X-rays to materialize as pairs in the gap. The replenished charges partially screen 
the electric field, which is self-consistently solved together with the distribution functions 
of particles and gamma-rays. If no current is injected at neither of the boundaries of 
the accelerator, the gap is located around the so-called null surface, where the local 
Goldreich- Julian charge density vanishes. However, we first find that the gap position 
shifts outwards (or inwards) when particles are injected at the inner (or outer) boundary. 
Applying the theory to the Crab pulsar, we demonstrate that the pulsed TeV fiux does 
not exceed the observational upper limit for moderate infrared photon density and that 
the gap should be located near to or outside of the null surface so that the observed 
spectrum of pulsed GeV fiuxes may be emitted via curvature process. 

Subject headings: gamma-rays: observation - gamma-rays: theory - magnetic field - 
pulsars: individual (Crab) - X-rays: observation 



1. Introduction 

The EGRET experiment on the Compton Gamma 
Ray Observatory has detected pulsed signals from 
seven rotation-powered pulsars (e.g., for Crab, Nolan 
et al. 1993, Fierro et al. 1998). The modulation of 
the 7-ray light curves at GeV energies testifies to the 
production of 7-ray radiation in the pulsar magneto- 
spheres either at the polar cap (Harding, Tademaru, 
& Esposito 1978; Daugherty & Harding 1982, 1996; 
Sturner, Dermer, & Michel 1995; Shibata, Miyazaki, 
& Takahara 1998), or at the vacuum gaps in the outer 
magnetosphere (Cheng, Ho, & Ruderman 1986a, b, 
hereafter CHR; Chiang & Romani 1992, 1994; Ro- 
mani and Yadigaroglu 1995; Romani 1996; Zhang & 
Cheng 1997, heafter ZC97). Effective 7-ray produc- 
tion in a pulsar magnetosphere may be extended to 
the very high energy (VHE) region above 100 GeV as 
well; however, the predictions of fluxes by the current 
models of 7-ray pulsars are not sufhciently conclu- 
sive (e.g., Cheng 1994). Whether or not the spectra 
of 7-ray pulsars continue up to the VHE region is a 
question which remains one of the interesting issues 
of high-energy astrophysics. 

In the VHE region, positive detections of radiation 
at a high confidence level have been reported from 
the direction of the Crab pulsar (Nel et al. 1993). 
However, as for pulsed TeV radiation, only the upper 
limits have been, as a rule, obtained. If the VHE emis- 
sion originates the pulsar magnetosphere, a significant 
fraction of them can be expected to show pulsation. 
Therefore, the lack of pulsed TeV emissions provides 
a severe constraint on the modeling of particle accel- 
eration zones in a pulsar magnetosphere. 

In fact, in CHR picture, the magnetosphere should 
be optically thick for pair-production in order to re- 
duce the TeV flux to an unobserved level by absorp- 
tion. This in turn requires very high luminosities of 
infrared photons. However, the required IR fluxes 
are generally orders of magnitude larger than the ob- 
served values (Usov 1994). We are therefore moti- 
vated by the need to contrive an outer-gap model 
which produces less TeV emission with a moderate 
infrared luminosity. 

High-energy emission from a pulsar magnetosphere, 
in fact, crucially depends on the acceleration electric 
fleld, along the magnetic fleld lines. It was Hi- 
rotani and Shibata (1999a, b,c; hereafter Papers I, II, 
III), and Hirotani (2000a,b,c; hereafter Papers IV, V, 
VI) who flrst solved the spatial distribution of to- 



gether with particle and 7-ray distribution functions. 
They demonstrated that a stationary gap is formed 
around the null surface at which the local Goldreich- 
Julian charge density, 

vanishes, where Bz is the component of the magnetic 
field along the rotation axis, il the angular frequency 
of the neutron star, and c the speed of light. Equation 
(|l|) is valid unless the gap is located close to the light 
cylinder, of which distance from the rotation axis is 
given by wj^c = c/fl. In this letter, we develop the 
method presented in Paper VI, by investigating the 
case when particles flow into the gap from the inner 
or the outer boundaries. 

In the next two sections, we describe the physical 
processes of pair production cascade and the resultant 
7-ray emission. We then apply the theory to the Crab 
pulsar and present the expected 7-ray spectra in § 4. 
In the final section, we compare the results with ZC97. 

2. Basic Equations and Boundary Conditions 

Let us first consider the Poission equation for the 
electrostatic potential, $. Neglecting relativistic ef- 
fects, and assuming that typical transfield thickness 
of the gap, Dj_, is greater than or comparable with 
the longitudinal gap width, W, we can reduce the 
Poisson equation into the form (Paper VI) 



ds 



— $ = 47re(iV+ - 7V_ - pcj/e) 



(2) 



where iV+ and iV_ designate the positronic and elec- 
tronic densities, respectively, e the magnitude of the 
charge on an electron, and s the length along the last- 
open fieldline. 

It is convenient to non-dimcnsionalize the length 
scales by a typical Debey scale length c/ujp, where 
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Be represents the magnetic field strength at the gap 
center. The dimensionless coordinate variable then 
becomes 

e ^ (^p/c)s. (4) 

By using such dimensionless quantities, we can rewrite 
the Poisson equation into 



(5) 
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where = e$(s)/(moC^); the particle densities per 
unit flux tube are defined by 
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~n B{s) ' 



(7) 



We evaluate B^/B at each point along the last-open 
field line, by using the Newtonian dipole field. 

We next consider the continuity equations for the 
particles. Assuming that both electrostatic and curvature- 
radiation-reaction forces cancel out each other, we ob- 
tain the following continuity equations 



±B^ 

ds\B 



1 



(8) 



where G±(s,e^) are the distribution functions of 7- 
ray photons having momentum ±moCe^ along the 
poloidal field line. Since the electric field is assumed 
to be positive in the gap, e^'s (or e^'s) migrate out- 
wards (or inwards). The pair production redistribu- 
tion functions, •q±^ are defined as 

'7p±(e7) (l^Mc)— / dex^7^crp(e7,ex,Mc), (9) 

where tip is the pair-production cross section and 
cos~^/ic refers to the collision angle between the 7- 
rays and the X-rays (see Paper VI for more details 
about eq. [^). 

Let us introduce the following dimensionless 7-ray 
densities in the dimensionless energy interval between 
/3j_i and Pi: 

9±i0^7T7r I de^G±{s,e,). (10) 



VlBc 



In this letter, we set /3o = lO'^, which corresponds to 
the lowest 7-ray energy, 51.1 MeV. We divide the 7- 
ray spectra into 9 energy bins and put fii = 10^'^, 
/32 = 103, /33 = 103-5, pi = 104, /35 = 104-5, pe = 
104-75, = iqS. Ps = 10^-^5^ and Pg = lO^-^. 
We can now rewrite equation (^) into 
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d^ 



i=l 



where rj^j. are evaluated at the central energy in each 
bin. 



A combination of equations ( |1 1| ) gives the current 
conservation law. 



jtot = "-+(0 + "-(C) = constant for ^. 



(12) 



When Jtot — 1-0, the current density per unit flux 
tube equals the Goldreich-Julian value, ^/{2tt). 

Unlike the charged particles, 7-rays do not prop- 
agate along the local magnetic field lines. How- 
ever, to avoid complications, we simply assume that 
the outwardly (or inwardly) propagating 7-rays di- 
late (or constrict) at the same rate with the magnetic 
field. This assumption gives a good estimate when 
W ^ njLc holds. We then obtain (Paper VI) 

^9±i0 = ^ (InS) T V:.±'9±' ± ^c^"±(0, (13) 
where i = 1, 2, • • •, m (to = 9) and 



Vc 



ds 



K5{t)dt, (14) 



where -f^s/s refers to the modified Bessel function of 
5/3 order. 

Equating the electric force e\d^/dx\ and the radi- 
ation reaction force, we obtain the saturated Lorentz 
factor at each point as 



Psat — 



3Rc' 
~27 



dx 



1 



1/4 



(15) 



we compute the curvature radius Rc at a point for 
a Newtonian dipole magnetic field. Since the maxi- 
mum of \d^/dx\ and the potential drop are roughly 
proportional to W'^ and W^, respectively (Paper V), 
the particles become unsaturated for very small W. 
To avoid an overestimation in such cases, we compute 
P by 

' ~ (16) 



P 



where (^(^2) represents the maximum attainable Lorentz 
factor. 

2.1. Boundary Conditions 

To solve the differential equations (||), (||), (pT|), 
and (|l|), we must impose boundary conditions. At 
the inner (starward) boundary (^ — d), we impose 
(Paper VI) 



S||(6) = 0, ^(6)=0, 



(17) 
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g+\^i)=0 (* = 1,2,...,9). 



(18) 



Since positrons may flow into the gap at ^ = as 
a part of the global current pattern in the magneto- 
sphere, we denote the positronic current per unit flux 
tube at ^ = ^1 as 



which yields (eq. [|2|) 



(19) 



(20) 



At the Oliver boundary (^ = ^2), we impose 
i?ll(6) = 0, 3-^(6) = (i = l,2,...,9), (21) 

^~(6)-j2. (22) 

The current density created in the gap per unit flux 
tube can be expressed as 



Jgap — Jtot — Ji — J2- 



(23) 



We adopt jgap, ji, and j2 as the free parameters. 

We have totally 24 boundary conditions (p^)-(|22 
for 22 unknown functions n±, g±'^ (z = 1, 2, • • 

•,9). Thus two extra boundary conditions must be 
compensated by making the positions of the bound- 
aries ^1 and ^2 be free. The two free boundaries ap- 
pear because E\\ = is imposed at both the bound- 
aries and because jgap is externally imposed. In other 
words, the gap boundaries (^1 and ^2) shift, if ji 
and/or j2 varies. 

3. TeV Spectra 

For simplicity, assume that the IR field arc homo- 
geneous and isotropic within the radius gjlc- Inter- 
polating the pulsed fluxes in radio, near IR, and opti- 
cal bands from the Crab pulsar (Moffett and Hankins 
1996; Percival et al. 1993; Eikenberry et al. 1997), 
we obtain (Paper V) 



'^^^^ = 1.5 X W^<f ^ 



IR 



V^LC 



(24) 



where tmracC^ refers to the IR photon energy, and 
eiR.min < e < eiR^max- Wc adopt eiR,min = 10"'' and 
eiR,min — 10^^; the results do not depend on these 
cut-off energies very much. 

If an electron or a positron is migrating with 
Lorentz factor F ^ 1 in an isotropic photon field. 



it upscatters the soft photons to produce the follow- 
ing number spectrum of 7-rays (Blumenthal & Gould 
1970): 



dtde-y 



3 c dNiji dem 

-0-rp- 



4 F2 deiR eiR 



2qlTLq+{l + 2q){l~q) + 



iQq)\l-q) 



2(1 + Qg) 



,(25) 



where Q = 4eiRF, q = / Q{T — e^); ctt is the Thom- 
son cross section, and the energy of the upscattered 
photons in rricC^ unit. Substituting equation (p^), in- 
tegrating d'^N/dtde^ over cir, and multiplying the 7- 
ray energy (e^TOoC^) and the electron number (iVe) in 
the gap, we obtain the flux density of the upscattered, 
TeV photons as a function of e^. 

We finally consider the extrinsic absorption of the 
TeV photons outside of the gap. For a homogeneous 
and isotropic IR field, the optical depth becomes 



r(e^) 



dN: 



IR 



rfeiR 



fp(eiR, £7, ^c)deiR, 



(26) 

where the path length is assumed to be ci7lc/2. 

For the Crab pulsar, r ~ 5 holds in TeV energies 
(fig. U). Therefore, the observed TeV flux reduces to 
about 1% of the intrinsic flux. 

4. Application to the Crab Pulsar 

In this section, we apply the theory to the Crab 
pulsar. The rotational frequency and the magnetic 
moment are fl = 188.1rad and /i = 3.38 x 
lO^OC cm^ 

4.1. Electric Field Structure 

HEAO 1 observations revealed that the X-ray spec- 
trum in the primary pulse phase is expressed by 



dN, 



pi 



A^pl^x (^min £x ^max). 



(27) 



with a = -1.81 and Npi = 5.3 x 10^''d'^{ro/TULc)~^ 
(Knight 1982), where d refers to the distance in 
kpc. We adopt e,„i„ = 0.1keV/511keV and e,„ax = 
50keV/511keV. Substituting this power-law spec- 
trum into equation (^ , we can solve the Vlasov equa- 
tions by the method described in § 2. 

We consider four representative boundary condi- 
tions: We choose (ji,j2) = (0,0), (0.3,0), (0.6,0), 
and (0,0.3) as cases 1, 2, 3, and 4, respectively. That 
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is, for case 2 (or case 4), the positronic (or electronic) 
current density flowing into the gap per unit flux tube 
at the inner (or outer) boundary is 30% of the typical 
Goldreich- Julian value, VI/2-k. We fix jgap = 0.01 for 
all the four cases, because the solution forms a 'brim' 
to disappear (fig. 2 in Hirotani & Okamoto 1998) if 
Jgap exceeds a few percent. In what follows, we adopt 
45° as the magnetic inclination. 

The results of E\\ {£) for the four cases are presented 
in figure ^. The abscissa designates the distance along 
the last-open field line and covers the range from the 
neutron star surface (s = 0) to the position where the 
disance equals s = 1.2zuj^c = 1.91 x 10® m. 

The solid line (case 1) shows that the gap is located 
around the null surface. However, the gap shifts out- 
wards as ji increases, as the dashed (case 2) and dash- 
dotted (case 3) lines indicate. This result is consis- 
tent with what was predicted in Shibata and Hirotani 
(2000) analytically. 

On the other hand, when j2 increases, the gap 
shifts inwards and the potential drop, ^{32), re- 
duces significantly. For example, we obtain $(52) = 
7.1 X 10^2 V for case 4, whereas 1.7 x 10^^ V for 
case 2. The reasons are sixfold: • In a station- 
ary gap, the pair production optical depth, W/Xp, 
equals the ratio iV^(jgap/jtot), where Ap and re- 
fer to the pair production mean free path, and the 
number of 7-rays emitted by a single particle, respec- 
tively. • The increased X-ray density at small radii 
reduces Ap. • The ratio jgap/jtot decreases as j2 in- 
creases. • As a result, W decreases very rapidly with 
increasing j2. • Owing to the rapidly decreasing W, 
$(52) significantly decreases, although the local pcj, 
and hence dE\\/ds increases at small radii. • As W 
decreases, decreases to some extent; however, this 
effect is passive and cannot change the conclusion. 

4.2. Gamma-ray Spectra 

The GeV spectra are readily computed from g+*(C2) 
and 5_*(^i), while the TeV spectra are obtained by 
the method described in § ||. We present the 7-ray 
spectra for the four cases in figure 0, multiplying the 
cross sectional area of D±^ = (6VFp. If D±^ increase 
twice, both the GeV and TeV fiuxes increases four 
times. 

In GeV energies, the observational pulsed spec- 
trum is obtained by EGRET observations (filled cir- 
cles; Nolan et al. 1993), while in TeV energies, only 
the upper limits are obtained by Whipple observa- 
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Fig. 1. — Pair production optical depth for a TeV 
photon to be aborbed in the homogeneous, isotropic 
IR field in the Crab-pulsar magnetosphere. 
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Fig. 2. — Distribution of E\^{S). The solid, dashed, 
dash-dotted, and dotted lines correspond to the cases 
1, 2, 3, and 4, respectively (see text). 
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tions (open squares; Weekes et al. 1989; Reynolds et 
al. 1993; Goret et al. 1993; Lessard et al. 2000), 
and by Durham observations (open triangle; Dowth- 
waite et al. 1984). The filled circles denote the un- 
pulsed flux obtained by CANGAROO observations 
(Tanimori et al. 1998). 

The figure shows that the TeV fluxes are kept be- 
low the observational upper limits as a whole for ap- 
propriate GeV fluxes. Therefore, we can conclude 
that the problem of the excessive TeV flux does not 
arise for a reasonable IR density for the Crab pulsar. 

It is noteworthy that the GeV spectrum depends 
on ji and j2 significantly. In particular, in case 4 
(as the dotted lines show), the GeV emission signifi- 
cantly decreases and softens, because both the poten- 
tial drop and the maximum of E'n reduce as the gap 
shifts inwards. As a result, it becomes impossible to 
explain the EGRET flux around 10^4 Hz, if the gap 
is located well inside of the null surface. 

5. Discussion 

In summary, we have developed a one-dimensional 
model for an outer-gap accelerator in the magne- 
tosphere of a rotation-powered pulsar. When the 
electronic current flows into the gap from the outer 
boundary, the gap shifts inwards to emit very soft 
GeV emissions. Applying this method to the Crab 
pulsar, we find that the gap should be located near 
to or outside of the null surface, so that the observed 
spectrum of pulsed GeV fluxes may be emitted via 
curvature process. By virtue of the absorption by the 
dense IR field in the magnetosphere, the problem of 
excessive TeV emission does not arise. 

Let us briefiy compare the present method with 
ZC97, who considered that the gap width is limited 
by the surface X-rays due to the bombardment of the 
particles produced in the gap. The magnetospheric 
X-rays considered in this paper is much denser than 
the surface X-rays due to the bombardment. As a 
result, the localized gap in the present paper produces 
less intrinsic TeV flux compared with what would be 
obtained in ZC97 picture. 

For cases 1,2, and 3, the intrinsic TeV luminosity 
is comparable or less than the GeV one. Therefore, 
the Lorentz factors are limited primarily by curvature 
process (eq.|l6)). For case 4, however, the intrinsic 
TeV luminosity well exceeds the GeV one; therefore, 
the radiation-reaction forces are due to IC scatter- 
ings rather than the curvature process. In fact, we 



may expect a sufficient GeV flux via IC scatterings 
when the gap is located well inside of the null sur- 
face. This is because the dense X-ray field will limit 
the particle Lorentz factors small (Paper II), and be- 
cause the less-energetic particles scatter copious IR 
photons into lower 7-ray energies with large cross sec- 
tions {'^ ctt). There is room for further investigation 
on this issue. 

One of the authors (K. H.) wishes to express his 
gratitude to Drs. Y. Saito and A. K. Harding for valu- 
able advice. He also thanks the Astronomical Data 
Analysis Center of National Astronomical Observa- 
tory, Japan for the use of workstations. 
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Fig. 3. — Gamma-ray spectra from the Crab pulsar 
magnetosphere. The thick (or thin) lines represent 
the flux of outwardly (or inwardly) propagating 7- 
rays. The solid, dashed, dash-dotted, and dotted lines 
correspond to the same cases as in figure 2. 
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